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ABSTRACT 


Lualualei  Naval  Magazine  is  located  in  the  northeastern 
corner  of  the  Lualualei  Valley,  Oahu,  Hawaii.  Tt  has  an  area 
of  approximately  12  square  miles  (between  latitudes  ?l029, 
and  2 1°2V25"  north  and  between  longitudes  158o09'  and 
158°05,^5"  west),  with  an  elevation  of  about  250  feet  above 
sea  level.  Shear  cliffs  of  the  Waianae  Range  bound  the  area 
to  the  north,  east  and  south,  and  alluvial  plans  of  the 
Lualualei  Valley  extend  all  the  way  to  the  raciflc  Ocean  to 
the  west.  During  late  1971*,  studies  were  carried  out  to 
determine  the  geothermal  potential.  Anomalously  low  values 
of  electrical  resistivities  were  mapped  in  the  southwest 
part  of  the  area.  In  addition,  shallow  (1  meter)  tempera- 
ture measurements  showed  that  the  area  with  anomalous  resis- 
tivities was  characterized  by  a temperature  of  ?(>.30C,  about 
2°C  above  the  normal  temperature  in  the  surrounding  area. 
Several  water  wells  drilled  in  the  1^50' s .also  indicate  that 
subsurface  temperatures  in  the  Lualualei  Valley  are  unusually 
warm.  A shallow  hole,  drilled  deep  enough  to  penetrate  the 
basalt  below  the  alluvium,  would  provide  a more  conclusive 
evaluation  of  the  prospect. 


INTRODUCTION 


The  islands  of  Hawaii  are  shield-shaped  basaltic  domes; 
therefore,  thermal  water  on  these  islands  owes  its  heat  to 
volcanic  sources.  It  may  well  be  that  the  rocks  in  Hawaii 
are  too  permeable  to  permit  steam  to  accumulate  under  pres- 
sure at  levels  close  to  the  surface,  and  the  abundance  of 
cold  ground  water  makes  unlikely  the  occurrence  of  steam  at 
a temperature  appreciably  above  the  boiling  point  (Macdonald 
and  Abbott,  1970).  A borehole  was  drilled  at  the  summit  of 
Kilauea  Volcano  during  the  summer  of  1973  (Keller,  1971!). 

This  study  proved  that  commercial  quality  steam  could  be 
obtained  at  a practical  depth.  Undoubtedly,  a tremendous 
amount  of  volcanic  heat  exists  at  relatively  shallow  depths 
in  the  islands  of  Hawaii. 

The  Waianae  caldera  on  the  island  of  Oahu  is  a late 
Tertiary  volcanic  center  which  exhibited  resurgent  activity 
during  Quaternary  period.  The  study  area  is  located  in  the 
Lualualei  Valley  which  is  situated  to  the  west  of  the  Waianae 
Range  (see  Figs.  1 anu  2).  Despite  the  possibility  of  exist- 
ence of  a heat  source,  there  are  no  geothermal  surface  mani- 
festations such  as  geysers,  hot  springs,  fumaroles,  solfatras, 
etc.  Considering  the  conditions  given  above,  this  geophysical 
investigation  was  directed  to  study  the  area  for  a buried  geo- 
thermal system. 

A review  of  geophysical  and  geologic  work  done  previously 
by  other  groups  was  made.  Gravity,  magnetic  and  seismic 
information  helps  define  the  geologic  environment  of  the 


igure  1.  Map  of  Oahu  showing  location  of  the  study  area 
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’"igore  2.  Picture  of  the  Luaiualei  Valley  showing  location  of  the  study  area 
(after  Macdonald  and  Abbott,  1970 ). 
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study  area  and  its  relationship  to  the  volcanic  center  and 
the  rift  zones. 

Electrical  resistivity  methods  have  been  successfully 
used  to  locate  geothermal  reservoirs  in  different  parts  of 
the  world  for  a number  of  years  (Keller,  1970).  In  order  to 
cover  the  area  in  detail  and  in  a short  period  of  time  the 
dipole  mapping  technique  (Furgerson,  1970)  was  chosen  as  the 
main  exploration  technique.  This  technique  is  very  effec- 
tive in  mapping  lateral  resistivity  variation  and  gives  a 
general  idea  of  the  vertical  variations.  Use  of  one  dipole 
source  is  not  enough  to  determine  best  apparent  resistivities 
and  anisotropies  can  cause  false  anomalies.  In  order  to 
determine  the  possibility  of  observing  a false  anomaly,  two 
dipole  sources  at  about  right  angles  to  each  other  were  set 
up  in  the  same  area.  Data  from  both  of  these  sources  were 
treated  together  using  the  rotating-quadripole  method  to 
compute  the  best  apparent  resistivities.  The  purpose  was  to 
delineate  an  area  with  low  resistivities,  inasmuch  as  geo- 
thermal reservoirs  contain  water  and  steam  at  high  tempera- 
tures and  pressures  and  dissolved  salts;  they  have  higher 
conductivities  than  the  surrounding  rocks. 

During  the  surveys,  a relatively  conductive  zone  was 
discovered.  A shallow  temperature  study  was  carried  out  and 
the  results  correlated  well  with  the  rotating-quadripole  data, 
increasing  the  likelihood  of  there  being  a sealed  geothermal 
reservoir  present. 
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The  Hawaiian  Islands  are  a chain  of  basaltic  shield  vol-  j 

canoes  built  over  a fissure  1600  miles  long  in  the  ocean  ! 

no  o t 

floor  between  154°40'  and  171°75'W  longitude  and  l8°54'  to  j 

28"15'N  latitude.  Recent  gravity,  magnetic  and  seismic  stud-  I 

ies  convincingly  showed  tnat  the  lavas  which  built  the  5 

Hawaiian  islands  were  extruded  primarily  along  faults  oriented 

1 

either  east-west  and  associated  with  the  Molokai  fracture  j 

system  or  northwest-southeast  and  associated  with  the  trend  I 

f 

along  which  lie  the  Koolau  dike  complex  and  the  Musician  j 

Seamounts.  The  volcanic  pipes  may  have  formed  at  points  of  j 

intersection  of  rifts  of  the  two  fracture  systems.  Whether  j 

the  NW-SE  fracture  system  is  a strike-slip  fault,  a simple  \ 

tension  crack  or  a tear  along  the  crest  of  a fold  is  not  * 

known.  To  the  southwest  of  Honolulu  at  about  1700  feet  below  j 

sea  level,  some  fossils  and  shallow-water  corals  of  late  j 

} 

Miocene  age  were  discovered  (MacDonald  and  Abbott,  1970).  j 

Therefore,  the  Hawaiian  Islands  might  have  been  built  on  an  f 

older  ridge,  perhaps  contemporaneous  with  a chain  of  islands  ] 

that  existed  in  middle  and  early  Tertiary  time  further  north-  | 

■t 

$ 

west.  The  NW-SE  fracture  system  probably  has  been  existent  l 

it 

\ 

since  early  Tertiary  (Stearns,  1966).  ; 

i 

The  Hawaiian  Islands  were  built  by  extrusive  materials  l 

from  volcanic  centers  during  late  Tertiary  and  Quaternary.  '] 


The  weight  of  the  extruded  lavas  caused  the  ridge  to  sink  in 
order  to  reestablish  isostatic  equilibrium.  Since  the  depth 
to  Moho  is  about  15  km  on  the  Ridge  and  about  11  km  in  the 
normal  ocean  basin,  a thickening  of  the  crust  of  some  *}  km 
is  indicated  (Strange,  Woollard,  and  Rose,  196*0.  Along  with 
volcanism,  erosion  has  taken  place.  Sea  cliffs  and  drainage 
systems  have  been  developed.  The  islands  have  continued  to 
submerge.  The  older  islands  to  the  northwest  are  believed  to 
have  submerged  up  to  about  9000  feet  (Stearns,  1966).  Con- 
currently, a new  epoch  of  volcanism  began  and  secondary  out- 
breaks continued  into  Recent  time.  The  eruptions  occurred  on 
all  the  major  islands  except  Lanai  in  Quaternary  period.  Com- 
plex submergences  and  emergences  continued.  Because  of  the 
rapid  alterations  of  these  events  very  little  reef  formation 
took  place. 

Recently  drilled  wells  in  Lualualei  Valley  penetrated 
1200  feet  of  stream-deposited  alluvium  before  entering  lava 
rock.  Since  the  bottom  of  the  valley  is  1200  feet  below  sea 
level  and  the  valley  had  to  be  cut  above  sea  level,  it  was 
concluded  that  the  island  must  have  sunk  at  least  1200  feet 
since  the  alluvium  was  deposited  (MacDonald  and  Abbott,  1970). 

Oahu  is  a volcanic  doublet  with  an  area  of  6j4  square 
miles  and  has  four  major  geomorphic  provinces.  They  are: 

1)  Koolau  Range,  2)  Waianae  Range,  3)  Schofield  Plateau,  and 
4)  Coastal  Plain  (see  Pig.  3). 
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The  Koolau  Range  is  the  younger  of  the  two  ranges  on 
the  island  and  it  is  located  to  the  east.  It  was  built 
principally  by  eruptions  along  a northwest-trending  rift  zone. 
The  rocks  of  the  Koolau  Volcano  are  mainly  tholeiltic 
basalts  and  olivine  basalts  with  small  amounts  of  oceanite. 
They  have  been  divided  into  two  groups.  The  Kailua  Volcanic 
series  are  the  eroded  rocks  of  the  ancient  Koolau  caldera. 


These  rocks  are  altered  by  hydrothermal  action  due  to  steam  * >, 

rising  in  the  vent  area.  The  Koolau  volcanic  series  are  ’ \ 

| 1 

those  lavas  and  dikes  lying  outside  the  caldera  and  are  < \ 

i i 
| i 

altered  only  rarely  by  hydrothermal  action.  After  the  Koolau  j .. 

1 \ 

volcano  ceased  activity,  a period  of  erosion  and  deposition  \ 

started  lasting  about  2 million  years.  Then  volcanic  j \ 

2 *J 

activity  resumed  on  the  southeastern  end  of  the  Koolau  range.  j '■ 

About  30  vents  have  erupted  and  they  are  called  the  Honolulu  \ 

volcanic  series.  The  .ents  are  aligned  chiefly  along  ' 

NE-SW  fissures  and  the  lavas  include  nephelinites , \ 
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member.  The  upper  member  is  the  relatively  thin  eap  that 
appears  to  have  covered  the  entire  top  of  the  shield  late  in 
its  history.  It  is  largely  hawaiite  with  lesser  amounts  of 
alkalic  olivine  basalt.  See  Table  1 for  stratigraphic  rock 
units  and  their  ages  (Stearns,  1966).  Some  posterosional 
eruptions  occurred  on  the  Waianae  Range,  during  the  Pleisto- 
cene, near  the  ancient  caldera.  These  alkalic  olivine 
basalts  are  called  Kolekole  volcanics  and  they  are  probably 
correlative  with  the  secondary  eruptions  on  the  Koolau  Range. 

The  Schofield  Plateau  was  formed  by  the  lavas  from  the 
Koolau  Range  banking  against  the  already-eroded  slope  of  the 
Waianae  Volcano  to  form  the  gently  sloping  surface  of  the 
Schofield  Plateau.  An  erosional  unconformity  between  the  rocks 
of  the  two  volcanos  is  visible  along  Kaukonahua  Gulch  at  the 
eastern  foot  of  the  Waianae  Range  (MacDonald  and  Abbott,  1970). 

The  coastal  plain  lies  mostly  on  the  ponded  lavas  of  the 
Koolau  Volcano  north  and  south  of  the  Schofield  Plateau.  The 
plain  is  composed  mainly  of  marine  sediments  deposited  on 
the  lavas  when  the  sea  stood  higher  In  mid-Pleistocene  time. 

See  Figure  4 for  the  geologic  map  of  the  island  of  Oahu. 


i 
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Geology  of  the  Study  Area 

After  the  Waianae  Volcano  ceased  activity,  stream  erosion 
began  and  great  valleys  were  carved,  especially  on  the  south- 
west side  of  the  Waianae  Range  where  the  streams  were  older 
and  the  rocks  weaker.  Then  the  island  went  through  a complex 
cycle  of  submergences  and  emergences  (Stearns,  1966).  The 


? 
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Table  1.  Stratigraphic  rock  units  in  the  island  of  Oahu  (after  Stearns,  1966). 
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Figure  4.  Geologic  map  of  the  island  of  Oahu  (after  Stearns,  1946) 
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great  submergence  resulted  in  deep  drowning  of  the  great 
vall-ys  and  their  subsequent  sedimentation.  Re-emergence  of 
Oahu  exposed  the  great  valleys  to  erosion  once  more,  with 
their  nearly  flat  alluvial  plains.  The  area  is  located  at 
the  northeast  corner  of  one  of  these  valleys,  namely  the 
Lualualei  Valley. 

The  study  area  is  covered  with  alluvium.  The  exact 
thickness  of  alluvium  in  the  vicinity  of  the  anomaly  is  not 
known;  however,  it  is  believed  to  be  about  1200  feet  (prob- 
ably less  than  1200  feet),  since  drill  holes  penetrated  1000 
feet  of  alluvium  near  the  town  of  Waianae  and  1200  feet  in 
the  Lualualei  Valley  (MacDonald  and  Abbott,  1970).  It  is 
expected  to  get  thinner  towards  the  steep  walls  of  the  valley, 
to  the  east,  the  north  and  the  south.  Outcrops  in  the  east- 
ern and  the  northeastern  part  of  the  area  are  believed  to  be 
from  the  Lower  Member  of  the  Waianae  Volcanic  Series.  Kole- 
kole  volcanics  which  was  the  result  of  renewed  volcanism  in 
the  Pleistocene  time  is  located  three  miles  to  the  northeast 
of  che  anomaly  (Pig.  *0 . The  caldera  of  the  Waianae  Volcano 
was  located  in  the  area  immediately  west  of  Kolekole  Pass, 
and  extended  from  the  northern  side  of  Makaha  Valley  to  the 
head  of  Nanakuli  Valley  (Fig.  5).  This  puts  the  study  area 
close  to  the  center  of  the  caldera.  Near  the  area  where  the 
Lower  and  the  Middle  members  of  the  Waianae  Volcanic  Series 
are  exposed,  there  are  numerous  dikes.  Their  sizes  vary 
between  a few  inches  and  15  feet,  and  most  of  them  are  nearly 
vertical  with  a general  trend  parallel  to  the  rift  zones. 


ubk  *****  > 


20  !0  15800  50  40  I 

Figure  5.  Shaded  relief  map  of  the  island  of  Oahu,  showing  location  of 

the  Waianae  Caldera  (after  Stearns,  A92&) 
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There  are  no  faults  large  enough  to  Justify  individual 
descriptions  in  or  near  the  study  area.  Numerous  dikes  of 
different  sizes  suggests  lateral  displacement  and  Jointing 
of  the  host  rock  with  each  injection.  Puu  Kailio,  which  i3 
a hill  located  at  the  head  of  Lualualei  Valley,  was  a place 
of  intense  volcanic  activity;  therefore,  it  was  decided  that 
the  Kailio  syncline  was  merely  a sag  in  the  flows  produced 
by  local  withdrawal  of  support  (Stearns,  1966).  The  anomaly 
is  located  about  2.5  miles  south-southwest  of  this  syncline 
(Fig.  2). 
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SURVEY  TECHNIQUES 


Rotatlng-Quadripole  Survey 

Description  of  the  Method:  The  rotating-quadripole  tech- 

nique is  an  electrical  resistivity  method  which  makes  use  of 
two  dipole  sources  oriented  at  approximately  right  angles  to 
each  other.  In  fact,  conventional  dipole  mapping  procedures 
were  used  in  the  field  (Furgerson,  1970,  and  Jordan,  197*0; 
however,  the  data  were  treated  differently. 

A suitable  location  was  selected  to  set  up  a fixed  d-c 
current  source,  keeping  in  mind  the  boundaries  oi  the  study 
area.  Then,  two  holes  800  meters  apart,  were  dug  to  bury 
two  tin  sheets,  3 feet  by  *1  feet.  In  order  to  reduce  the 
contact  resistance,  salt  and  water  wsre  used.  These  two 
source  electrodes  were  connected  with  #12  Insulated,  copper 
wire.  To  generate  an  electric  field,  a 2.5  kw  generator  was 
used.  The  1-phase,  120v  a-c,  60  hz  output  of  the  generator 
was  passed  through  an  asymmetrically-timed  switching  circuit 
and  then  rectified  to  d-c  current.  The  asymmetric  square-wave 
output  of  the  rectifier  was  fed  into  the  source  wire.  The 
relative  polarity  of  the  voltage  drop  along  receiver  wire  was 
determined  by  observing  the  asymmetry  of  the  square-wave  output. 

The  field  measurements  were  made  a quarter  of  a mile  apart 
from  each  other  along  the  roads  and  near  the  ammunition  maga- 
zines. At  every  station,  two  components  (preferably  ortho- 
gonal) of  the  earth  potential  were  measured.  To  do  this, 
two  non-polarizing  receiver  electrodes  were  placed  30  meters 
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apart  and  connected  to  a d-c  amplifier  with  #12  insulated 
copper  wire.  The  output  of  the  amplifier  was  recorded  with 
a potentiometric  chart  recorder.  A Brunton  Compass  was  used 
to  determine  azimuths  of  both  component:..  Measurements  were 
made  at  103  stations  around  the  first  source.  Then  the 
second  source  was  set  up  at  a 96°  angle  to  the  first  source. 
Measurements  were  repeated  occupying  the  same  stations. 

Two  voltage  components  for  each  source  were  divided  by 
the  receiver  length  (30  m)  to  obtain  electric  field  components 
and  added  vectorially  to  determine  the  total  electric  fields 
for  both  sources.  If  both  sources  were  turned  on,  there  would 
be  only  one  total  electric  field  and  it  can  be  obtained  by 
orthogonal  addition  of  the  total  electric  fields  of  the  two 
sources.  The  fields  from  the  two  sources  can  also  be  added 
with  appropriate  weights  to  rotate  the  total  electric  field 
l80°.  A set  of  apparent  resistivities  for  every  station  was 
calculated  using  the  rotating  quadripole  equation  (see  Appendix 
A-l  for  the  equation,  and  Appendix  A-2  for  the  data  obtained 
in  the  field).  Out  of  40  resistivity  values  the  best  (the 
highest)  value  was  selected  for  every  station. 

Interpretation  of  Field  Data 

Best  Apparent  Resistivity:  The  dipole  mapping  technique 

has  been  used  extensively  to  determine  lateral  resistivity 
variations  due  to  conductive  bodies  at  depth,  such  as  geo- 
thermal reservoirs.  This  fast  d-c  resistivity  method  also 
reveals  information  about  the  geologic  structure  of  an  area 
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provided  a resistive  basement  is  present,  and  can  be  used  for 
depth  estimates.  However,  the  technique  is  not  problem  free. 

Presence  of  a fault-like  boundary  between  the  region  with 
moderate  resistivities  and  another  region  with  high  resis- 
tivities will  cause  "false"  anomalies  to  be  observed  at 
specific  locations  along  the  fault-like  boundary  (Pig.  6) 

(Grose  and  Keller,  197*0. 

To  eliminate  the  observed  false  anomalies  and  to  obtain 
resistivities  as  close  as  possible  to  the  true  resistivities 
multiple  coverage  must  be  provided.  This  can  be  done  by 

* 

randomly  located  multiple  sources  or  by  the  rotating-quadripole  j 

method.  In  the  rotating  quadripole  method,  two  dipole  sources  j 

i 

at  right  angles  to  each  other  are  used.  j 

•i 

T 

Out  of  40  resistivity  values  obtained  by  rotating  the 

i 

total  electric  field,  the  highest  value  was  picked  for  every 

i 

station.  These  values  are  believed  r,o  be  the  closest  to  the 
true  resistivities;  therefore,  they  were  called  the  best 
apparent  resistivities.  In  almost  all  cases  the  best  resis- 
tivities are  higher  than  the  resistivities  obtained  from  the 
original  dipole  source. 

Figure  7 shows  the  best  apparent  resistivity  map  of  the 
study  area.  The  7b’  ohm-m  contours  and  high  resistivities 
indicate  the  presence  of  a resistive  basement  at  depth.  The 
resistive  basement  is  believed  to  be  the  volcanic  plug  of  the 
Waianae  Volcano.  A closed  35  ohm-m  contour  near  the  sources 
indicates  the  presence  of  a thin  conductive  surface  layer. 
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Figure 


7. 


Best' apparent  resistivity  map  of  the  Lualuclei  Naval  Magazine  — rota 
quadripole  arroy 
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The  35  and  30  ohm-m  contours  to  the  southwest  of  the  sources 

fr 

might  be  caused  by  a conductive  body  at  depth. 

A histogram  showing  the  distribution  of  best  resistiv- 
ity values  is  presented  in  Pig.  8.  A deviation  from  a normal 
distribution  occurs  between  28  and  4l.i|  ohm-m,  and  these 
anomalous  values  were  measured  mainly  in  the  conductive  zone  j 

located  to  the  southwest  of  the  two  sources  (Magazines  11  and  j 


12). 

To  construct  a scatter  diagram,  hest  resistivities  were 
plotted  against  the  distance  to  the  nearest  of  the  four  source 
electrodes.  This  distance  is  thought  to  be  the  depth  of 
current  penetration.  The  resulting  scatter  diagram  (Pig.  9) 
suggests  the  presence  of  layering  in  the  area,  because  the 
points  tend  to  cluster  along  a ^5°  angle  after  1 km  distance. 


Temperature  Survey 

Geothermal  exploration  should  not  be  limited  to  the 
vicinity  of  areas  of  heat  escape  in  the  form  of  geysers, 
fumaroles,  sulfatras  and  hot  springs.  A reservoir  of  moder- 
ate size  with  an  upper  surface  at  2 km  depth  would  approxi- 
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mately  double  the  normal  geothermal  gradient  over  an  area 
of  a few  square  kilometers  (Banwell,  1970).  Therefore, 
detection  of  a sealed  reservoir  with  no  convection  to  the 
surface  is  not  difficult.  Surface  temperature  and  heat  flow 
measurements  at  1 to  2 m depths  have  been  successful  in 
mapping  high  temperature  zones  and  in  evaluating  the  ho#t  water 
distribution  below  the  surface.  Such  observations  provide  a 
rapid  and  direct  way  of  making  an  estimate  of  the  size  and 


Apparent  Resistivity 


Figure  9.  Best  apparent  resistivity  scatter  diagram 
rotating-quadripole  array 
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energy  potential  of  a system.  When  making  or  interpreting 
a shallow  temperature  survey,  hydrological  conditions  in 
the  area  including  rate  and  direction  of  ground  water  move- 
ment have  to  be  considered.  Because  movement  of  ground 
water  can  carry  away  the  conductive  heat  flow  from  a 
thermal  anomaly,  it  can  displace  the  surface  temperature 
patterns.  Instrumental  errors  may  result  from  using  thermis- 


tor probes  for  temperature  measurements;  the  overall  error 

-O 


has  been  estimated  to  be  between  0.01  and  0.05  C for  meas- 
urements at  a depth  of  2 m (Dedkova,  et  al,  1970).  Seasonal 
and  weather  changes  and  variations  in  surface  cover  are 
believed  to  have  random  effects  on  shallow  temperatures. 

Description  of  the  Method:  A shallow  (1  meter  deep) 

temperature  study  was  carried  out  in  the  Lualualei  Naval 
Magazine  in  order  to  determine  the  shallow  temperature  dis- 
tribution and  locate  areas  with  anomalous  temperature  gra- 
dients. A simple  mechanical  device  called  a dongeyknocker 
made  of  mild  steel  was  used  to  punch  holes  in  the  ground. 

A precision  thermistor  teflon  probe  was  placed  at  the  bottoms 
of  these  holes  using  a wooden  rod.  The  resistance  of  the 
probe  was  measured  with  a digital  multimeter.  These 
resistances  were  converted  to  temperatures  (°C)  using  a con- 
version table.  The  measurements  were  spaced  closely  (800 
feet  apart)  over  the  resistivity  anomaly  to  provide  detail 
in  that  area.  A wider  spacing  of  2000  feet  was  used  for  the 
measurements  over  the  resistive  zone.  A total  of  70  measure- 
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ments  was  made  including  7 to  the  west  of  the  Pence  Road 
where  there  is  no  resistivity  information. 

Interpretation  of  the  Field  Data 

Ground  Temperature:  Measured  temperature  values  repre- 

sent ground  temperatures  and  Pig.  10  shows  the  ground  tem- 
perature (at  1 meter  depth)  pattern  of  the  area.  Before 
making  an  attempt  to  interpret  the  data,  some  important  fac- 
tors should  be  considered.  Surface  cover  in  the  area  is  allu- 
vium and  major  variations  likely  to  cause  false  anomalies  are 
unlikely.  Since  there  is  a difference  of  2°C  in  temperatures 
between  the  cold  and  the  warm  zones,  instrumental  errors 
(.05°C  or  less)  could  not  have  changed  the  temperatures 
enough  to  distort  the  temperature  pattern.  Where  ground  water 
is  confined  between  dikes,  it  Is  under  static  conditions. 

This  Is  true  for  the  basalts  with  numerous  dikes  below  the 
alluvial  fill;  however,  the  alluvial  section  itself  is 
believed  to  contain  waver  under  dynamic  conditions. 

It  was  discovered  that  the  area  with  anomalous  resis- 
tivities (less  than  30  ohm-m)  had,  in  general,  a temperature 
of  26.3°C.  Since  the  resistive  zones  had  lower  ground  tem- 
peratures (as  low  as  21J.3°C),  it  was  concluded  that  the  low 
resistivities  to  the  southwest  of  the  dipoles  were  caused  by 
high  temperatures  at  depth,  possibly  a sealed  geothermal 
reservoir.  The  other  temperature  anomalies  in  the  area 
(areas  with  temperatures  higher  than  26°C)  are  believed  to 
be  surfleial  features  caused  either  by  variations  in  the 
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Figure  10*  Ground  temperature  map  of  the  Lucluaiei  Naval  Magazine,  based  on 
measurements  mode  at  1 m depth 
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surface  cover  or  by  the  movement  of  the  ground  water  in 
alluvium  or  both,  and  they  are  not  supported  by  low  resis- 
tivities except  where  there  is  no  resistivity  information. 

Self-Potential  Survey 

Geothermal  reservoirs  contain  ground  water  of  differing 
temperature  and  chemical  composition.  Thermal  gradients  in 
pore  water  electrolytes  and  contact  potentials  between  bodies 
of  ground  water  of  differing  temperature  and  chemical  com- 
position may  give  rise  to  measurable  electrical  anomalies. 

In  the  course  of  resistivity  surveys  of  thermal  areas,  high 
natural  electrical  potentials  have  been  commonly  observed 
(Banwell,  1970).  These  natural  potentials  may  form  a diag- 
nostic pattern  over  the  area.  Such  a survey  would  be  valu- 
able since  it  would  have  some  depth  penetration  and  point  by 
point  observations  would  be  possible  with  small  electrode 
spacings. 

Description  of  the  Method:  A self-potential  survey  was 

carried  out  In  order  to  determine  if  the  area  with  low 
resistivities  and  high  temperatures  had  anomalous  self  poten- 
tials. The  measurements  were  made  with  a digital  voltmeter 
using  two  non-polarizing  electrodes  with  an  electrode  spacing 
of  200  m.  A starting  point  with  an  assigned  potential  of 
zero  was  picked  in  an  area  with  high  resistivities  (Pig.  11). 
Potential  drops  were  measured  along  a continuous  line,  across 
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the  area  of  interest,  ending  in  an  area  with  high  resistiv- 
ities. A loop  was  made  around  the  anomaly  to  distribute  the 
total  error  among  individual  measurements . 

Interpretation  of  Field  Data 

Self  Potential:  After  error  distribution  (2  mv/200  m), 

potential  drops  with  respect  to  the  initial  point  with  zero 
potential  were  calculated.  Contoured  self-potential  values 
are  shown  in  Pig.  11.  Self-potential  over  the  area  decreases, 
in  general,  to  the  north  and  there  is  no  indication  of  an 
anomalous  pattern  in  the  area  with  low  resistivities  and  high 
temperatures. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Does  the  area  where  the  Lualualel  Naval  Magazine  is 
located  have  a geothermal  reservoir?  If  so,  does  it  contain 
high  enthalpy  water  at  an  economical  depth?  These  were  the 
questions  under  consideration  during  the  course  of  this  study. 

The  area  is  located  near  a center  of  intense  volcanism 
(Puu  Kailio)  inside  the  volcanic  vent  zone  of  the  Waianae 
Caldera.  Even  though  the  Waianae  Range  was  built  in  late 
Tertiary,  the  volcanism  continued  during  the  Quaternary  period. 

Young  volcanism  in  the  area  suggests  the  presence  of  a strong 
and  active  heat  source.  All  geophysical  work  done  in  the  area 
indicates  the  presence  of  a dense  body  of  rock  extending  from 
shallow  depths  (800  .n  from  the  surface)  to  a depth  of  5.5  km 
under  the  Waianae  Caldera.  It  is  probable  that  olivine  and 
perhaps  other  heavy  crystals  that  have  lagged  behind  In 
basaltic  magma  rising  toward  the  surface,  and  have  accumu- 
lated to  form  a plug  of  peridotite  that  resembles  mantle  rock 
in  its  physical  properties.  The  drill  hole  at  the  summit  of 

Kilauea  penetrated  rocks  with  low  permeability  and  densities  j 

3 

around  3.2  gr/cc.  The  volcanic  plug  under  the  Waianae  Cal- 
dera is  believed  to  be  made  up  of  rocks  similar  to  the  rocks 
found  under  the  Kilauea  Volcano.  Low  permeability  probably 
minimized  the  heat  escape  from  the  heat  source;  thus,  maintain- 
ing its  heat  content  over  the  years.  The  geological  condi- 
tions in  the  area  indicate  the  presence  of  a heat  source. 
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The  rocks  located  above  the  volcanic  plug  are  believed 
to  be  highly  fractured  due  to  a high  number  of  dikes  of  dif- 
ferent sizes.  Since  these  highly  permeable  basalts  are 
located  above  the  possible  heat  source,  and  since  the  area 
has  abundant  groundwater  due  to  high  rainfall,  the  section 
between  the  volcanic  plug  and  the  alluvial  fill  has  the  prop- 
erties to  behave  as  the  reservoir  rock  for  a geothermal  system. 
The  area  where  the  anomaly  is  located  may  have  up  to  1000  feet 
of  alluvium  on  the  surface.  The  alluvium  (Pleistocene)  is 
believed  to  be  partially  consolidated  and  inasmuch  as  it  is 
located  above  the  permeable  basalts,  and  there  is  a possibil- 
ity of  hydrothermal  alteration,  the  alluvium  could  behave  like 
a caprock  because  of  its  reduced  permeability.  Presence  of  a 
zone  with  low  permeability  over  the  reservoir  minimizes  the 
heat  escape  to  the  surface,  keeping  most  of  the  heat  in  the 
reservoir,  thus  increasing  the  temperature  of  the  groundwater. 
From  the  discussion  above,  it  can  be  concluded  that  the  geo- 
logic environment  is  suitable  for  presence  of  a geothermal 
system. 

To  locate  the  area  with  anomalous  geothermal  gradient, 
a rotating-quadripole  survey  was  carried  out . Anomalous  low 
resistivities  were  discovered  to  the  southwest  of  the  dipole 
sources.  A shallow  temperature  study  which  was  done  after- 
wards confirmed  that  the  low  resistivities  were  caused  by 
high  temperatures  in  the  area.  Because  there  is  no  surface 
activity,  chemical  analyses  to  determine  the  reservoir  tem- 
perature were  not  available.  Temperature  measurements  made 
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at  a meter  depth  could  not  be  used  to  determine  the  true  geo- 
thermal gradient . 

As  a result  of  this  study,  it  is  concluded  that  the  area 
where  magazines  11  and  12  are  located  may  have  a sealed 
geothermal  reservoir  at  depth.  However,  the  vertical  resis- 
tivity distribution  arid  the  temperature  gradient  in  the  area 
are  not  known.  A test  hole  which  is  deep  enough  to  penetrate 
the  basalts  below  the  alluvium  should  be  drilled  to  determine 
the  temperature  gradient  and  also  to  obtain  valuable  informa- 
tion about  the  rocks  in  the  area.  Perhaps  a number  of  initial 
test  holes  with  intermediate  depths  (20  m)  could  be  drilled  to 
determine  the  hottest  part  of  the  anomalous  area. 


as 


34 


Take  two  dipole  sources  at  approximately  right  angles  to 
each  other.  Even  though  it  is  desirable,  they  do  not  have  to 
have  a common  electrode. 


The  potential  expression  for  apparent  resistivity  assum- 
ing spherical  current  spreading  from  single  electrode  in  a 
homogeneous  earth  is  u * pl/2irB  (Van  Nostrand  and  Cook,  19 66). 
The  potential  at  point  P due  to  all  four  electrodes  can  be 
written. 


U2  « -pi;L/2Tr c (X-Xx ) 2+  Cy-y^^  5 2i1/2 


U2  - +pl1/2ir[(x-x2)2+(y-y2)231/2 


U3  * -pl2/2ff[(x-x3)2+(y-y3)2]1/2 


- +pl2/2x[(x-xi|)2+(y-y4)2]1/2 


Since  E * -VU,  E - ,-3U/3x  and  E * -3U/3y 

x y 


x and  y components  of  the  electric  field  at  point  P due  to  all 


four  electrodes: 


Eix  * P^Cx^^/auCCx-x^^Cy-yi)2]372 


Ely  “ Pli(y“yi)/2wC(x-x1)2+(y-y1)233/2 


E2x  “~PIi^x“x2^/2,,r^X”X2^2+^y“y2^233/2 


E2  *-pI1(y-y2)/2ir[(x-x2)  +(y-y2)  3 


2.3/2 


E3x  " pI2(x-x3)/2^C(x-x3)2+(y-y3)233^2 


E3y  “ PI2^“y3>/2,r^x-x3>2+(y"y35^3/2 


Enx  *"pI2  ^ x~x4 )/2fff ( x“x4 ) 2+ Cy~y n ) 2 3 3/2 


E/jy  a,-pI2(y-yi1)/2iT[(x-x24)2+(y-y^)233/2 


The  total  electric  fields  in  x and  y directions  are: 


*Tx  2tt  1 


1l  ^x"xi^ 


I1  (x-x2) 


[(x-x^+fy-y^^372  [(x-x2)2+(y-y2)23‘ 


I2  (x-X3) 


IgCx-x^) 


[(x-x3)2+(y-y3)233/2  [(x-xlj)2+(y-yij)233/2 
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[(x-x1)‘!+(y-y1)‘:];,/‘:  C(x-x2)2+(y-y2)z3- 


i2  <y-y3)  i2  (y“y4) 

+ C<x-x3)5+(y-y3>¥7S  ‘ [(x-x4)5+(y-y,)233/5} 


The  amplitude  of  the  total  electric  field  at  point  P is: 


p p 1/2 

E_._  » (Em  + Em  ) 

TOT  v Tx  Ty ' 


Let: 


Ri  " £<x-xi>2  + (y-yx)2]172 


R2  - C(x-x2)2  + (y-y2)231/2 


r3  - [(x-x^)2  + (y-y^)2]172 


Ru  - [(x-x,)2  + (y-y,.)2]172 


Substitute  four  equations  above  in  the  equation  for  E„ 


'’TOT  2 it 


CjCx-Xi)  Ij^Cx-Xg)  I2(x-x3)  Igfx-x^)  2 
* + . — ] 


rIl(y_yl)  Il(’  X2)  I2(y“y3)  I2(y“y4).2.1/2 

+ [ r + * J } 

RJ  q 3 p 3 ^ 

^ **Jj 


-r\:  t'fyt.'fe  '£'**'T>v  ' ' ; 

&?Li  L •,  i-i.fi.r  k y oJl&A  L.jM't&S  ,'A 


I,  (x-X,  ) I,  (X-X,)  I,(X~X,)  I,(x-x„) 


iUX  clt 


21^ (x-x^) (x-x2)  2I1I2(x-x1) (x-x,) 

vV  + ^ 

-2I1I2(x-:  (x-Xjj ) -2I1I2(x-x,)  (x-x,) 

R2 

2 

21,  I„(x-x,)  (x-x,. ) 21,  (x-x_)  (x-X|, ) 


p 3r  3 
R2  r4 


h3\ 


i12(y -yx>2  i12(y-y2)2  i2z(y-y,) 

+ — t + r + 


I22(y-Y1|)  2I12(y-y1)(y-y2) 


V V 


2I1I2(y-y1)(y-y3)  2I1I2(y-y1 ) (y-y^) 

4»  Tirmi*.  ,i.  -n-  1 ...  ■*  r ‘i  i.*i  ' ^L-m 


RlV 


2I1I2(y-y2)(y-y3)  2I1l2(y-y 


R23r33 


r23f 
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2i2  (y-y3)(y-ylj)  x/2 
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I 2‘  I 2 

ETOT  * fir  {ZT  C(x-Xi)2+(y-y1)2]  + t(x-x2)2+(y-y2)23 

R1  R2 


I c X 

+ C(x-x3)2+(y-y3)2]  + -g-  [ (x-x4)2+(y-y4)2] 

r3  r4 


■3 — ^ C(x-x1)(x-x2)  + (y-y1)(y-y2)] 


R1  R2 


2I1X2 


R1  R3' 


C(x-x1)(x-x3)  + (y~y1>(y-y3)] 


2IlX2 


R1  1 V 


[(x-x^Cx-Xj,)  + (y-y1)(y-yi|)] 


2I1I2 


VV 


t(x-x2)(x-x3)  + (y-y2)(y-y3)] 


2IiI2 

+ — 5 — j [(x-x2)(x-x4)  + (y-y2)(y-y4)] 
R2  r4 


2I2  1/2 

— ^ ^ C(x-x3)(x-x4)  + (y-y3)(y-y4)]> 
r3  r4 


Let:  C * I1/I2  and  substitute  this  equation  in  the  last  equa- 


tion for  ETQrp  above. 


-jttia, s .W.4& t5aaii_ .ft- anu.  * ,.x iuw.,  a . ».^.w  « ’*s$trL 


* 


P1?  1 1 

ETOT  - w-  'rr  ♦ rr  + pr  + pr 

n2  n3  n4 


w 


C(x-x1)(x-x2)  + (y-y1)(y-y2)] 


+ — ^ ^ [(x-x1)(x-x3)  + (y-y1)(y-y3)] 
R1  R3 


3 — j C(x-x1)(x-x1<)  + (y-y1)(y-y1{)] 


R1  R4 


— ^ ^ [(x-x2)(x-x3)  + (y-y2)(y-y3)] 

r2  r3 


+ — 3 — 3 [(x-x2)(x-x4)  + (y-y2)(y-y4)3 
R2  Ri< 

2 1/2 

3 — I C(x-x3)(x-X2|)  + (y-y3)(y-y1{)]> 

R3  Ri^ 


Call  that  part  of  the  equation  in  parentheses  K. 


ry  _ £ l r . _ 

xOT  2it  P I„K 


The  expression  for  apparent  resistivity  is 


where  K Is  the  geometric  factor.  To  find  the  apparent  resis- 
tivity expressions  for  dipoles  one  and  two  individually,  let 
Ip  = 0 and  1^  = 0 respectively  In  the  composite  equation  for 
the  rotating-quadrlpole  array  above. 


Afo  fry*,  ,, 
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Let: 


h 9 0 


Pi 


“TOT,  - 2T  <7?  + rr  " 7§T3  CC^Kx-x^^y-y^Cy-y^]} 

**2  *^X  *^2 

p-^i 

ETOT1  “ 2 iT  K1 
2irETOT1 

p * — y-nr—  : Apparent  resistivity  equation  for  the  first 
1lnl 


1/2 


source. 


Let: 


h ' 0 


pi- 


{— ij|-  + — ^ — j[(x-x3)(x-x1|)+(y-y3)(y-y1|)]} 


E « __2  r 1 
TOT_  2 x o-1  p -»  p jn 

d «3  Kg  K|j 


pi. 


1/2 


•a  mi  Y 

aT0T2  2tr  *2 

2it£ipQfpp 

p » — __ — ; Apparent  resistivity  equation  for  the  second 
x2*2 


source. 
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Station  number 


Distance  from  measuring  point  to  the  origin  (the  inter- 
section point  of  the  two  sources) (km) 


BEAR: 


CUR1 


CUR2 


Rmax 


Rmin 


Azimuth  of  the  line  between  the  origin  and  the  receiver 
station  (degrees) 

Azimuth  of  the  first  receiver  line  (source  #1) (degrees) 
Azimuth  of  second  receiver  line  (source  #1) (degrees) 
Azimuth  of  first  receiver  line  (source  #2) (degrees) 
Azimuth  of  second  receiver  line  (source  §2) (degrees) 


Voltage  measured  along  first  receiver  line  (source  #l)(mv) 
Voltage  measured  along  second  receiver  line  (source  #l)(ipv) 
Voltage  measured  along  first  receiver  line  (source  #2)  (mv) 
Voltage  measured  along  second  receiver  line  (source  #2)(mv) 
Total  step  current  on  source  #1  (amperes) 

Total  step  current  on  source  #2  (amperes) 


Length  of  receiver  line  (meters) 
Resistivity  due  to  source  #1  (ohm-m) 
Resistivity  due  to  source  #2  (ohm-m) 
Maximum  resistivity  (ohm-meters) 


Minimum  resistivity  (ohm-meters) 


Ratio  of  maximum  to  minimum  resistivity 
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APPENDIX  A- 3 

TEMPERATURE  STUDY  FIELD  DATA 
N:  Station  number 
R:  Resistance  (ohms) 

T:  Temperature  (°C) 


APPENDIX  A- 4 


SELF-POTENTIAL  STUDY  FIELD  DATA 


N:  Station  number 

V:  Potential  drop  along  a 200  m line  (millivolts) 

Vc:  Corrected  value  of  V (after  error  distribution) 
(millivolts) 

VN:  Potential  drop  between  a station  and  the  initial 
point  with  assumed  zero  potential  (millivolts) 


' ^ ^ »•  '-i I VW * < ' r’ U A - vfe. W J»  - «*  _‘,w 


-23.9 
-33.4 
- 6.2 

- 7.9 
-27.7 

24.5 

-54.5 

- 2.9 
-13.2 
-24.6 

8.0 

-19.6 

- 4.7 

- 9.9 

0.2 

-12.8 

- 8.1 

- 7.7 
-11.3 
-31.0 

5.2 

- 8.4 

-32.2 

- 5.8 
-16.1 

- 4.8 

- 1.5 

- 7.8 
-15.9 
-24.0 
-40.0 
-30.8 
-15.1 


-21.9 

-31.4 

- 4.2 

- 5.9 
-25.7 

26.5 

-52.7 

- 0.9 
-11.2 
-22.6 

10.0 

-17.6 

- 2.7 

- 7.9 
2.2 

-10.8 

- 6.1 

- 5.7 

- 9.3 
-29.0 

7.2 

- 6.4 
-30.2 

- 3.8 
-14.1 

- 2.8 

0.5 

- 5.8 
-13.9 
-22.0 
-38.0 
-28.8 
-13.1 


-21.9 

-53.3 

-57.5 

-63.4 

-89.1 

-62.6 

-115.3 

-116.2 

-127.4 

-150.0 

-140.0 

-157.6 

-160.3 

-168.2 

-166.0 

-176.8 

-182.9 

-188.9 

-197.9 

-226.9 

-219.7 

-226.1 

-256.3 

-260.1 

-274.2 

-277 

-276.5 

-282.3 

-292.2 

-318.2 

-356.2 

-385.0 

-398.1 


N 

V 

vo 

VN 

34 

3.8 

5.8 

-392.3 

35 

-20.2 

-18.2 

-410.5 

36 

-1.4 

0.6 

-225.5 

37 

10.8 

12.8 

-214.7 

38 

- 2.7 

- 0.7 

-215.4 

39 

10.8 

12.8 

-202.6 

40 

12.7 

14.7 

-187.9 

41 

13.4 

15.4 

-172.5 

42 

8.1 

10.1 

-162.4 

43 

18.5 

20.5 

-141.9 

44 

8.3 

10.3 

-131.6 

45 

11.2 

13.2 

-118.4 

46 

30.8 

32.8 

- 75.6 

47 

-30.0 

-28.0 

-103.6 

48 

- 1.0 

1.0 

-102.6 

49 

43.2 

45.2 

- 67.4 

50 

- 2.7 

- 0.7 

-258.0 

51 

27.3 

29.3 

-228.7 

52 

-12.0 

-10.0 

-238.7 

53 

6.8 

8.8 

-229.9 

54 

1.6 

3.6 

-226.3 

55 

0.8 

2.8 

-222.5 

56 

- 0.8 

1.2 

-221.3 

57 

10.6 

12.6 

-208.7 

58 

-10.2 

- 8.2 

-216.9 

59 

- 8.9 

- 6.9 

-223.8 

60 

0.3 

2.3 

-221.5 

61 

- 3.4 

- 1.4 

-222.9 

62 

0.1 

2.1 

-220.8 

63 

-19.0 

-17.0 

-237.8 

64 

- 4.0 

- 2.0 

-239.8 

65 

6.0 

8.0 

-231.8 

50 


0 

0 


REFERENCES 


Banwell,  C.  J.,  1970,  Geophysical  techniques  in  geothermal 
exploration:  UN  Symposium  on  the  Development  and 

Utilization  of  Geothermal  Resources,  Pisa,  1970,  sec.  IV. 

Dedkova,  D.,  Halousek,  J.,  Krcmar,  B.,  and  Prihoda.-  K. , 1970, 
Geothermal  prospection  in  shallow  holes  and  its  limita- 
tion: UN  Symposium  on  the  Development  and  Utilization 

of  Geothermal  Resources,  Pisa,  1970,  IV/5. 

Furgerson,  R.  B. , 1970,  A controlled-source  telluric  current 
technique  and  its  application  to  structural  investiga- 
tions: Thesis  1313*  Colo.  School  of  Mines,  Golden, 

Colo. 


I 


.‘i 


0 

0 

0 

[j 

G 


0 

0 


Furgerson,  R.  B. , and  Keller,  G.  V.,  1974,  Computed  dipole 
resistivity  effects  for  an  earth  model  with  vertical 
and  lateral  contrasts  in  resistivity:  Report  No. 

NR081275,  Office  of  Naval  Research,  Washington,  D.C. 

Grose,  L.  T.,  and  Keller,  G.  V.,  1974,  The  Colorado  School  of 
Mines  Nevada  Geothermal  Study  - Report  of  Progress  for 
the  period  May  1,  1974  to  July  31,  1974:  Colo.  School 

of  Mines,  Golden,  Colo. 

Jordan,  J.  M. , 1974,  Geothermal  investigations  in  the  San 

Luis  Valley,  South-Central  Colorado:  Thesis  1478,  Colo. 

School  of  Mines,  Golden,  Colo. 

Keller,  G.  V.,  1974,  Drilling  at  the  summit  of  Kilauea  Volcano: 
Prepared  for  National  Science  Foundation,  Colo.  School 
of  Mines,  Golden,  Colo, 

1966,  Dipole  method  for  deep  resistivity  studies: 

Geophysics,  v.  31,  no.  6,  p.  1088-1104. 

1970,  Induction  methods  in  prospecting  for  hot 

water:  Geothermics,  spec,  issue  2,  p.  318-332. 

Macdonald,  G.  A.,  and  Abbott,  A.  T.,  1970,  Volcanoes  in  the 
sea,  the  geology  of  Hawaii:  Honolulu  Hawaii,  University 

of  Hawaii  Press. 

Stearns,  H.  T. , 1966,  Geology  of  the  State  of  Hawaii:  Palo 

Alto,  Calif.,  Pacific  Book3. 

Strange,  W.  E,,  Woollard,  G,  P.,  and  Rose,  J.  C.,  1965,  An 
analysis  of  the  gravity  field  over  the  Hawaiian  Islands 
in  terms  of  crustal  structure:  Pacific  Sci.,  v.  19, 

no.  3,  p.  381-389. 

Van  Nostrand,  R.  G.,  and  Cook,  K.  L. , 1 966.  Interpretation  of 
resistivity  data:  U.S.  Geol.  Survey  *rof.  Paper  499 • 


J 


APPENDIX  B 


r 


U 

U 

u 

[J 


EXPLORATION  ON  ADAK  ISLAND, 


ALASKA 


by 


David  L. 


Butler  and  George  V.  Keller 


u 

Li 

U 


( 


i 

u 


i 

; 


i 


Appendix  B accompanies  a report  entitled  "Geothermal  Energy 
in  the  Pacific  Region"  by  L.  T.  Grose  and  G.  V.  Keller, 

May,  1975.  This  project  has  been  supported  by  the  Off i..  a 
of  Naval  Research  Contract  Number  N00014-71-A-0430-0004. 


APPENDIX  B - EXPLORATION  ON  ADAK  ISLAND,  ALASKA 
David  L.  Butler  and  George  V.  Keller 


Introduction 

As  pointed  out  in  the  body  of  this  report,  the  Aleutian  Islands 
comprise  an  area  of  active  to  recently  active  volcanism  which  should 
be  favorable  for  the  occurrence  of  geothermal  energy.  Adak  Naval 
Base  is  an  excellent  candidate  for  the  use  of  geothermal  energy, 
inasmuch  as  Adak  Island  is  located  toward  the  western  end  of  the 
active  portion  of  the  Aleutian  Islands.  Because  of  this,  some  ex- 
ploration was  carried  out  on  Adak  Island  in  an  attempt  to  evaluate 
the  geothermal  potential.  Seismicity  and  resistivity  surveys  were 
carried  out  on  Mount  Adaydak,  on  northern  Adak,  near  the  Naval  Base 
(see  Figure  1) . 

The  area  of  interest,  the  northern  mountainous  part  of  Adak, 
is  the  remnants  of  three  historically  inactive  volcanoes.  From  east 
to  west,  the  volcanoes  are  Mount  Adagdak,  Andrew  Bay  Volcano,  and 
Mount  Moffet.  Although  no  extensive  faulting  has  been  mapped  in 
this  area,  some  evidence  of  faulting  is  present  on  the  north  slope 
of  Mount  Adagdak.  The  lithology  of  this  area  includes  volcanic  and 
intrusive  rocks  of  Paleozoic  age.  Tertiary  to  Quaternary  volcanic 
rocks  and  Quaternary  alluvium.  The  alluvium  unit  consists  largely 
of  glacial  drift  and  other  unconsolidated  materials,  including  vol- 
canic ash. 


Seismicity  Surveys 
o 

More  than  1000  km  near  Adak  Island  was  surveyed  from  Oct.  22 
to  November  1,  1974,  for  microearthquakes  to  aid  in  the  evaluation 
of  the  geothermal  potential  of  the  area.  The  survey  was  carried 
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out  by  Microgeophysics,  Inc.,  under  contract  to  the  Colorado  School 
of  Mines,  The  size  of  the  1000  km  area  is  based  on  a detection 
threshold  of  magnitude  0 or  less  with  at  least  one  station  record- 
ing an  event  within  the  area.  The  objective  was  to  detect  and 
locate  microearthquakes  and  thereby  map  tectonically  active  struc- 
tures. 

Figure  2 is  an  historical  seismicity  map  of  the  area 
around  the  Adak  Island.  The  Aleutian  Arc  is  an  active  seismic 
province,  however,  a large  number  of  the  events  shown  on  Figure 
2 occur  on  or  near  the  subduction  or  Benioff  Zone  of  the 
Aleutian  Island  arc-trench  system.  The  subduction  zone  is  the 
postulated  tectonic  feature  occurring  where  one  crustal  plate 
is  being  forced  under  another.  In  the  Aleutian  Islands,  the 
North  Pacific  Plate  is  being  moved  in  a northwesterly  direction 
under  the  American  Plate. 

The  subduction  zone  strikes  east-west  south  of  Adak  and 
dips  to  the  north  to  a depth  of  100  to  150km  beneath  the  island. 
The  seismic  activity  associated  with  the  subduction  zone  occurs 
at  too  great  a depth  to  be  of  interest  in  this  survey  but  the 
tectonic  activity  at  depth  is  a manifestation  of  the  regional 
stress.  This  same  stress  pattern  may  influence  surface  faulting 
and  the  shallow  microseismicity  of  interest  in  this  survey. 

The  specific  project  area  contains  no  historical  epicenters 
although  several  epicenters  are  northeast  of  the  area  of  interest. 
In  early  March,  1974,  13  small  earthquakes  were  detected  by  the 
NOAA  Seismological  Observatory  on  Adak  and  located  near  Andrew 
Bay  and  Mount  Adagdak  (Mr.  D.  Glover,  personnel  communication, 
1974).  These  events  are  of  small  magnitude  (m<3)  and  are  not 
plotted  on  the  historical  seismicity  map  (Figure  2) . 
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During  the  course  of  this  survey,  assistance  and  data  were 
supplied  by  Mr.  D.  Glover,  observatory  director  of  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA) , Seismological 
Observatory,  Adak,  Alaska.  The  data  supplied  includes  arrival 
times,  S-P  times  and  magnitudes  of  several  local  earthquakes. 

His  enthusiastic  assistance  is  gratefully  acknowledged. 

The  next  section  of  this  paper  outlines  the  instrumentation 
and  operational  methods  employed  in  the  field  work.  The  obser- 
vations and  analysis  are  given  in  a section  and  are  followed  by 
an  interpretation  of  the  results.  Recommendations  are  listed 
in  the  last  section  of  the  body  of  the  report.  The  appendix  is 
a listing  of  times,  locations  and  magnitudes  of  the  local  earth- 
quakes detected  during  this  survey. 

Instrumentation  and  Operational  Summary 

Seven  Sprengnether  Instrument  Co.  MEQ-800-B  portable 
seismic  systems  were  used  for  this  survey.  Each  system  con- 
sists of  a Mark  Products  model  LC-4,  1-hz  natural-frequency 
vertical  seismometer,  gain-stable  amplifier,  integral  timing 
system,  and  smoked  paper  recording  with  0.025mm  stylus  width 
and  120mm/min  recording  speed.  The  frequency  characteristics 
of  the  instrument  are  summarized  in  Figure  3 (note  that  both 
the  velocity  response  and  the  displacement  response  are  plotted; 
displacement  response  at  a particular  frequency,  f,  is  obtained 
by  multiplying  the  velocity  response  by  2irf ) . Gain  changes  are 
by  6db  steps  from  the  typical  operating  level  of  +96db  plotted 
in  the  figure. 
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Clocks  were  synchronized  daily  with  a master  clock,  which 
was  synchronized  with  WWV.  Clock  drifts  between  synchronizations 
were  below  expected  record  reading  errors,  therefore  no  correct- 
ions were  necessary.  Records  were  read  to  +0.03  sec  for  P 
arrivals  and  +0.10  sec  for  S-P  times.  Amplitudes,  peak  to  peak, 
were  read  to  the  nearest  millimeter,  and  durations  to  the 
nearest  0.5  sec. 

Despite  the  poor  working  conditions,  considerable  effort 
was  made  to  locate  stations  on  hardrock  outcrops  (crystalline 
exposures  or  well  compacted  sediments) . All  stations  were 
operated  at  the  gain  limit  allowed  by  ambient  background  noise. 
Station  locations  of  the  array  operated  for  ten  days  are  illus- 
trated in  Figure  4 and  listed  in  Table  1.  Stations  1 through  9 
were  operated  by  MicroGeophysics  Corporation?  zl-z5  and  AD8  by 


NOAA.  Details  of  the  operation  of  stations  1 through  9 are  shown 
in  Table  2.  The  NOAA  stations, which  are  operated  continuously, 
are  part  of  the  permanent  tsunami-warning  network  of  NOAA.  The 
operation  schedule  (Table  2)  shows  the  gains  and  the  time  periods 
of  individual  station  occupations. 
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STATION  COORDINATES 
TABLE  1 


Station 

X*  in  km 

Y*  in  km 

Z*  in  km 

1 

1.33 

7.60 

0.00 

2 

6.62 

4.95 

0.00 

3 

5.38 

6.03 

0.03 

4 

5.15 

11.75 

0.09 

5 

4.80 

11.40 

0.09 

6 

4.15 

8.42 

0.03 

7 

7.31 

8.58 

0.05 

8* 

0.00 

0.00 

0.09 

9 

3.70 

4.60 

0.03 

z2 

1.70 

0.35 

0.06 

z3 

0.05 

-1.90 

0.03 

z4 

-1.60 

0.15 

0.09 

z5 

0.15 

0.05 

0.09 

AD8 

5.60 

11.00 

0.24 

* The  origin  is  located  at  station  8 (latitude  51.88°N,  longitude 
176.69°W).  Positive  X and  Y are  east  and  north  respectively. 

Z is  the  station  height  above  sea  level  (altitude) . 
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Observations  and  Analysis 

Events  were  regarded  as  seismic  in  origin  if  they  appeared 
on  two  or  more  stations  with  time  differences  corresponding  to 
seismic  velocities  or  if  they  were  similar  in  appearance  to 
other  larger  events  which  were  well  recorded.  Seismic  events 
were  considered  local  if  they  had  S-P  times  of  less  than  four 
seconds. 

The  S-P  time  is  a characteristic  indicator  of  distance  to 
the  epicenter  of  an  event.  The  S-P  time  is  the  difference 
between  the  arrive!  time  of  the  S or  shear  wave  and  the  P or 
compressional  wave.  The  two  body  waves  propagate  by  different 
mechanisms  and  at  different  velocities  dependent  on  the  para- 
meters of  the  transporting  medium.  The  S-P  time  is  thus  a 
function  of  the  distance  traveled. 

An  example  of  a local  event  near  Adak  is  shown  in  Figure  5. 
Regional  and  teleseismic  (distant)  events  with  S-P  times  greater 
than  four  seconds,  were  considered  outside  the  scope  of  this 
survey  and  therefore  no  attempt  was  made  to  locate  them.  The 
four-second  S-P  time  cut-off  for  local  events  is  an  arbitrary 
limit  chosen  by  the  interpreter. 

Local  events  timed  on  four  or  more  stations  were  located 
using  a generalized  inverse  computer  program.  This  program 
assumes  a velocity  model  and  least-squares  fits  calculated 
travel  times  to  the  observed  arrival  times.  The  velocity 
structure  in  the  project  area  can  be  estimated  from  a layered 
crustal  velocity  model  of  the  Adak  region  developed  by  the  USGS 
(Engdahl,  1974) . Figure  6 illustrates  the  USGS  model  and  two 
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other  velocity  structures,  a half-space  model  and  a linear 
velocity -increase -with-depth  model.  The  half-space  model  was 
used  to  test  the  picked  arrival  times  and  the  sensitivity  of 
the  locations  to  velocity  changes.  The  constant  velocity  model 
obtained  good  fits  with  no  anomalous  station  residuals.  The 
velocities  obtained  with  the  half-space  model  agree  with  the 
averaged  velocities  of  the  layered  model.  The  model  that 
produced  the  best  fits  and  lowest  residuals  in  the  project  area 
was  the  linear  velocity -increase -with -depth  model.  This  model 
is  in  general  agreement  with  the  layered  model  and  it  produced 
station  residuals  on  the  order  of  the  picking  error  (+.03  sec). 

Figure  7 shows  all  located  events  and  illustrates  that  the 
local  seismicity  is  confined  to  two  areas.  Events  from  these 
two  areas  were  identified  by  using  S-P  times  and  signature 
similarities  and  were  classified  as  "type  A"  or  "type  B"  events. 
Figure  8 shows  a plot  of  the  cumulative  number  of  events  record- 
ed during  the  survey  of  each  type.  Type  A events  were  recorded 
at  a rate  of  four  events  per  day,  type  B events  at  a rate  of 
three  per  day  when  there  was  activity. 

Type  B events  were  located  30  to  40km  southeast  of  Adak. 

Due  to  the  distance  from  the  array,  the  epicenter  location 
precision  is  about  +4km.  Because  these  events  are  well  outside 


the  area  of  interest,  no  further  analysis  was  performed. 

Twenty  six  events  denoted  as  type  A were  located  in  the 
project  area.  The  precision  of  location  for  these  events  is 
about  +lkm  in  plan  and  +2km  in  depth.  The  events  occur  within 
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Andrew  Bay,  west  and  northwest  of  Mount  Adagdak.  The  located 
type  A events  are  shown  in  Figure  9.  The  least-squares-f it 
fault  plane  of  the  event  location  is  shown  in  Figure  10.  The 
least-squares-f it  solution  is  arrived  at  in  the  following 
manner.  After  the  surface  fault  trace  is  plotted  in  plan  view, 
an  auxilary  or  cross  section  view  is  taken  perpendicular  to 
the  strike  of  the  fault  trace  (to  show  the  true  dip  of  the 
fault  plane) . The  depth  perspective  of  the  fault  is  found  by 
plotting  the  depths  of  epicenters  (hypocenters)  in  the  auxilary 
view.  After  allowing  for  uncertainties  in  the  depth  parameter 
from  the  velocity  model  (+2km) , a fault  plane  can  be  fitted  to 
the  cross  section.  The  fault  plane  strikes  N60°E  and  dips  NW 
at  70°.  The  dotted  line  is  the  expected  precision  of  the  fault 
plane  solution. 

The  direction  of  movement  on  the  fault  plane  is  obtained 
by  the  use  of  a first  motion  study.  This  solution  is  shown  on 
Figure  11.  The  first  motion,  compression  for  vertical  motion 
up  and  dilatation  for  vertical  motion  down,  is  plotted  on  the 
upper  hemisphere  of  a sterographic  projection.  A fault  plane 
is  then  fitted  to  the  data  by  dividing  the  plot  into  four 
quadrants.  The  shaded  region  on  Figure  11  represents  compression 
and  the  unmarked  region  represents  the  dilatation  quadrants.  The 
first-motion  method  produces  two  orthogonal  solutions.  If  you 
chose  the  fault  plane  solution  which  strikes  N70°E  with  a Nw 
dip  of  75°  the  relative  motion  of  the  fault  is  right-lateral, 
strike-slip  with  a small  compones  c of  thrust.  The  fault  plane 


ADAK  FIRST  MOTION 


o DILATATION 
• COMPRESSION 


■'V. 

M*  • 

1 u 

» »v 
"\  • 
' \ ! 
\'  ■ 


'\\  *^>'o 


,•  ••  \\<s 

• **  Us 

..••:**  V\' 


mffi 0% 

i noPc 


AMPLITUDE  VS  DISTANCE 


solutions  from  the  two  methods  agree 
well  within  reasonable  limits  of 
accuracy. 

To  estimate  the  level  of  seis- 
micity some  attention  must  be  paid 
to  the  size  or  magnitude  of  an 

CO 

earthquake.  In  this  survey  an  0 

/ 0) 

amplitude  versus  distance  plot  (see 

< 

Fiqure  12)  was  constructed  for  two  ]> 

ch  - 2 

local  events  (October  30  at  6“j/ 
and  at  8h15m) . The  magnitude  of 
these  events  was  determined  by  using 
curves  established  for  microearthquakes 
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bv  Brune  and  Allen  (1967) . Magnitudes  from  two  other  larger 
events  (October  28,  llh04m  and  October  30  at  8h03m)  were  supplied 
by  the  NOAA  Adak  Observatory.  These  four  events  were  then  used 
to  determine  a magnitude  versus  duration  curve  (see  Figure  13) . 

To  correct  for  local  geologic  effects  all  durations  were  measured 
at  station  3. 

From  the  magnitude  versus  duration  curve,  all  detected 
events  were  assigned  a magnitude.  The  resulting  cumulative 
recurrence  curve  (a  plot  of  the  log  of  the  number  of  events 
of  a given  magnitude  versus  the  magnitude)  is  shown  in  Figure  14. 
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"he  slope  3f  this  linear  relationship  i often  called  the  b-slope. 

A b-slope  of  -0.8  is  shown  with  the  data  for  reference. 

Figure  15  shows  the  relationship  between  Poisson’s  ratio  and 
thv  velocity  ratio  . The  often  assumed  value  of  Poisson's  ratio 
of  ”.35  is  shown,  but  an  increase  of  up  to  0.1  often  occurs  when 
...  tcck  is  fractured.  The  velocity  ratio  is  defined  as  the  ratio 
of  the  compressional  wave  velocity  (Vp)  divided  by  the  shear  wave 
''olocity  (Vs).  Poisson's  ratio  is  defined  as  the  ratio  of  the 
transverse  contraction  to  the  longitudinal  extension  of  a rod 
runiected  to  an  axial  load.  Poisson's  ratio  is  also  a qualitative 
measure  of  the  amount  of  fracturing  present  in  a rock  mass, 
therefore  a map  indicating  variations  of  Poisson's  ratio  or  the 
\ ..  uOO’..ly  ratio  may  indicate  areas  or  volumes  with  anomalous 
fracturing  and  resultant  high  permeability. 

Figure  16  is  a plot  of  arrival  time  versus  S-P  times  for 
: rant  events  recorded  by  this  survey.  This  relationship,  often 
called  the  Wadati  diagram  method,  has  been  discussed  by  others 
•.Jemyenov,  A.N.,  1969;  Nersesov,  et.al.,  1971,‘  Kisslinger, 
n'-'dahl,  1973).  These  eight  events  showed  good  S-wave  breaks 
■ nd  i . e well  recorded  across  the  net.  Figure  16  gives  two  values 
; i each  slope.  The  A parameter  is  the  slope  of  the  linear  re- 
. a i. 1 .^i'lship  which  indirectly  gives  the  velocity  ratio  and  a gives 
' .v.  urresponding  Poisson's  ratio  from  Figure  15. 

Figure  17  is  a plot  of  the  velocity  ratio  with  time. 

Ja.  iu-iens  from  a constant  velocity  ratio  have  been  used  to 
_ r'.g  , -sren)onitory  phenomena  for  large  earthquakes  with  some 


THIS  FIGURE  SHOWS  THE  RELATIONSHIP  BETWEEN 
POISSON'S  RATIO  AND  VELOCITY  RATIO.  THE 
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success  (Semyenov,  1969;  Nersesov,  1971;  Kisslinger,Engdahl,1974) 
The  short  time- sampling  period  of  this  survey  precludes  any  pre- 
diction from  this  data  but  the  spatial  distribution  of  Poisson's 
ratio  can  be  interpreted. 
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To  locate  an  anomalously  fractured  zone,  a map  of  the 
observed  Poisson's  ratio  (plotted  at  the  epicenter)  was  made 
(see  Figure  18) . Note  that  the  form  lines  indicate  trends  in 
the  data  and  are  not  contour  lines.  The  Poisson's  ratio  was 
assumed  to  be  more  heavily  influenced  by  the  source  region  and 
mechanism  than  by  the  travel  path,  therefore  the  value  was 
assigned  to  the  epicenter.  Assignment  of  the  value  to  the  travel 
path  would  complicate  the  display,  but  the  same  conclusion  would 
be  derived  from  the  alternate  display.  Figure  18  produces  a 
consistent  picture  of  Poisson's  ratio  increasing  to  the  west  and 
with  depth  along  the  fault. 
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The  anomalous  distribution  of  Poisson's  ratio  and  its 
spatial  relation  to  the  geologic  features  can  be  explained  by 
a theoretical  model.  This  model  is  calculated  by  assuming  a 
Volterra-type  dislocation  for  the  fault  mechanism.  For  this 
model,  the  volume  strain  or  dilation  can  be  calculated  around 
the  fault  given  the  spatial  ceometry  of  the  fault  plane.  The 
Volterra-type  dislocation  model  assumes  that  a homogeneous  half- 
space is  broken  along  a distinct  plane  and,  after  displacement, 
is  welded  back  together.  The  half-space,  which  was  initially 
stress  free,  now  is  subjected  to  a regional  stress  (manifested 
by  the  volume  strain,  which  can  be  contoured)  created  by  the 
dislocation  at  the  "fault  plane".  The  dilation  or  volume  strain 
around  such  a dislocation  model  has  been  calculated  by  Yeatts 
(1975) . Figure  19  shows  a plan  and  cross  section  view  of  the 
dilation  around  a right-lateral,  strike-slip  fault  dislocation. 
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The  shaded  areas  represent  areas  under  dilitation  and  unshaded 
areas  are  under  compression.  The  region  of  most  severe  volume 
expansion  can  be  expected  to  exhibit  microfractures,  demonstrated 
by  an  increase  in  Poisson's  ratio,  leading  to  increased  perm- 
eability. 

Resistivity  Survey 

A resistivity  survey  was  carried  out  in  the  Mount  Adagdak 
peninsula  area  during  the  period  October  25  to  November  5,  1975,  by 
a party  under  the  supervision  of  Dr.  Paul  Donaldson  of  the  Colorado 
School  of  Mines.  The  rotating  dipole  technique,  as  described  by 
Tasci  in  Appendix  A of  this  report,  was  used.  Operations  in  the 
field  on  Adak  in  late  Fall  are  frought  with  difficulties,  largely 
induced  by  the  cold,  wet,  and  windy  weather.  Relatively  few  meas- 
urements were  made  before  the  survey  had  to  be  terminated,  and  the 
resistivity  survey  must  be  considered  to  be  only  preliminary  in 
nature. 

In  the  rotating  dipole  method,  field  data  are  obtained  by  mak- 
ing measurements  of  electric  field  intensity  around  pair  of  bipole 
sources  located  at  a single  site,  but  with  different  orientations. 
For  the  Adak  survey,  a pair  of  sources,  each  700  meters  in  length, 
v/as  placed  on  the  southwest  slope  of  Mount  Adagdak , as  shown  on  the 
map  in  Figure  21.  Currents  of  a few  amperes  were  used  to  excite 
these  two  source  bipoles.  This  is  less  than  normally  used  in  dipole 
mapping  surveys,  but  only  a small  motor  generator  set  was  available 
for  the  survey,  and  grounding  resistances  were  surprisingly  high. 
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Measurements  were  made  at  the  receiver  sites  indicated  on 
Figure  21,  using  receiver  lines  100  meters  in  length.  Voltages 
were  recorded  on  paper-tape  recorders.  Values  of  apparent  resistiv- 
ity were  computed  for  each  source  independently,  and  for  various 
combinations  of  the  two  sources  to  yield  resistivity  values  as  a 
function  of  current  flow  directional  at  each  receiver  site.  The 
maximum  value  of  apparent  resistivity  is  plotted  and  contoured  on 
Figure  21,  because  this  value  is  relatively  immune  to  problems 

caused  by  false  anomalies. 

The  roughly  elliptical  pattern  of  contours  shown  on  this  map, 
with  high  values  of  apparent  resistivity  being  observed  close  to 
the  source,  is  characteristic  of  an  earth  in  which  resistivity 
decreases  with  depth.  It  appears  that  at  distances  of  less  than 
1 kilometer  from  the  source,  current  has  not  penetrated  deep  enough 
to  be  controlled  by  the  properties  of  rocks  below  sea  level . At 
the  greater  distances,  beyond  1 kilometer,  apparent  resistivity 
values  decrease  to  less  than  20  ohm-meters,  or  values  which  are  com- 
patible with  the  existence  of  geothermal  fluids  at  depth.  However, 
much  more  detailed  coverage  of  the  area  will  be  required  before 
the  existence  of  any  potential  geothermal  reservoir  can  be  indicated. 


Interpretation 

A seismically  active  zone  exists  near  Andrew  Bay,  north 
and  west  of  Mount  Adagdak  on  Adak  Island,  Aleutians  Islands, 
Alaska.  This  zone  is  interpreted  as  a right-lateral,  strike- 
slip  fault  with  a small  component  of  thrust.  The  surface  trace 
of  this  fault  should  be  apparent  on  the  west  side  of  Mount 
Adagdak. 

A zone  of  relatively  high  Poisson's  ratio  mapped  in  this 
area  indicates  that  a fractured  or  dilated  region  exists  in  the 
area.  Consideration  of  the  microearthquake  station  locations 
and  the  travel  paths  from  the  recorded  events  indicates  the 
dilation  zone  probably  exists  south  of  and  on  the  west  end  of 
the  active  fault  zone. 

Reconciliation  of  the  theoretical  model  with  the  field 
observations  also  indicates  that  a dilated  zone  should  be  found 
north  of  and  on  the  eastern  end  of  the  active  fault  plane.  Due 
to  the  array  geometry  and  location  of  the  active  seismic  region, 
no  data  was  obtained  to  confirm  or  deny  the  existance  of  this 
second  anamolous  zone. 

Due  to  the  increased  permeability  caused  by  the  fracturing 
associated  with  the  mechanism  of  faulting  in  the  region,  the 
hypothesized  zones  of  volume  expansion  along  the  active  trace 
of  the  fault  are  interpreted  as  having  the  highest  potential  of 
producing  commercial  earth  steam.  These  areas  are  shown  in 
Figure  20.  The  interpretation  is  of  course,  subject  to 
clarification  when  reconciled  with  additional  geological  and 
geophysical  data  in  the  project  area. 
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